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q \ Abstract. The primary objective of the KATRIN experiment is to probe the absolute 

" ^ ' neutrino mass scale with a sensitivity of 200 meV (90% C.L.) by precision spectroscopy 

of tritium /3-decay. To achieve this, a low background of the order of 10 -2 cps in the 
region of the tritium /3-decay endpoint is required. Measurements with an electrostatic 
retarding spectrometer have revealed that electrons, arising from nuclear decays in the 
volume of the spectrometer, are stored over long time periods and thereby act as a 
major source of background exceeding this limit. In this paper we present a novel 
Q\ . active background reduction method based on stochastic heating of stored electrons 

\ by the well-known process of electron cyclotron resonance (ECR). A successful proof- 

of-principle of the ECR technique was demonstrated in test measurements at the 
KATRIN pre-spectrometer, yielding a large reduction of the background rate. In 
addition, we have carried out extensive Monte Carlo simulations to reveal the potential 
£SJ ■ of the ECR technique to remove all trapped electrons in a few ms with negligible 

loss of measurement time in the main spectrometer. This would allow the KATRIN 
experiment attaining its full physics potential. 
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The Karlsruhe Tritium Neutrino (KATRIN) experiment [I] is a next generation 
direct neutrino mass experiment, currently under construction at Tritium Laboratory 
Karlsruhe at the Karlsruhe Institute of Technology (KIT) Campus North site. KATRIN 
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is designed to measure directly the effective electron anti-neutrino mass m Ve , defined as 



where U e i denotes elements of the Pontecorvo-Maki-Nakagawa-Sakata leptonic mixing 
matrix and nij denote the neutrino mass eigenstates [2] . The KATRIN design sensitivity 
of 200 meV at 90% confidence level will allow to cover the quasi-degenerate neutrino 
mass pattern (where ni] » m 2 « m 3 ) and to investigate the role of relic neutrinos from 
the Big Bang in the evolution of large-scale structures in the universe [I]. 

1.1. The KATRIN setup 

The experiment analyzes the shape of the tritium /3-decay spectrum in a narrow region 
close to its endpoint at Eq = 18.6 keV. A non-zero neutrino mass reduces the maximum 
energy of the electron and changes the shape of the tritium /3-spectrum in the immediate 
vicinity of the endpoint. To reach the neutrino mass sensitivity a high energy resolution 
spectrometer, a high signal count rate as well as a low background rate are required. 

The 70 m long KATRIN setup (see figure [1]) combines a molecular windowless 
gaseous tritium source (WGTS) of highest stability and luminosity [3] with a high 
resolution electrostatic spectrometer for precision /3-spectroscopy. A magnetic guidance 
system directs the electrons created in the tritium source to the spectrometer section, 
consisting of a smaller pre-spectrometer providing the option to filter out low- 
energy electrons, and a larger main spectrometer for precision energy analysis. Both 
spectrometers work as electrostatic filters and transmit only those electrons that have 
sufficient kinetic energy to pass the retarding potential. The transmitted electrons are 
counted at a segmented silicon (Si-PIN) detector. By measuring the count rate for 
different filter voltages, the shape of the integral energy spectrum close to the tritium 
endpoint can be determined and m Ve can be deduced. 

To collimate the momenta of the /3-electrons created isotropically in the WGTS, 
the magnetic field drops by four orders of magnitude from the entrance/exit of 
the spectrometer to its center, the so called analyzing plane. The combination of 
electrostatic filtering with magnetic adiabatic collimation is called the MAC-E-Filter 
principle [HE] and allows for both large solid angle acceptance and high energy 
resolution. 

1.2. Background due to stored electrons 

The magnetic field configuration of a MAC-E-Filter generically creates a magnetic 
bottle for light charged particles being produced in the volume of the spectrometer. 
When propagating towards the entrance or exit parts of the spectrometer, they are 
magnetically mirrored by the increasing magnetic field (see figure [2]) [&H8] . Sources of 
stored electrons with energies up to several hundred keV stem from a-decays of the 
radon isotopes 219 ' 220 Rn and /3-decays of tritium molecules (HT, T2) occurring in the 
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Figure 1. Main components of the KATRIN experimental setup, a: Rear section, 
b: windowless gaseous tritium source, c: differential pumping section, d: cryogenic 
pumping section, e: pre-spectrometer, f: main spectrometer, g: focal plane detector. 
Below, the magnetic field and the electric potential along the beam axis are displayed. 
In both spectrometers the MAC-E filter principle is applied: As the electric potential 
is increased to U re t = —18.6 kV to filter the /3-electrons, the magnetic field drops from 
B c = 6 T to B c = 3 • 1(T 4 T. The magnetic flux of $ = 191 Tern 2 is conserved along 
the beam line. 

volume of the spectrometer. Stored electrons in the multi-keV-range mainly lose energy 
via ionization of residual gas molecules thereby generating secondary electrons at lower 
energies [9HlTj. Due to the excellent ultra-high vacuum (UHV) conditions [T2l[T3"] at 
the KATRIN spectrometers, the storage time of the keV-range electrons can last up to 
several hours. In this time period hundreds of secondary electrons are produced. 

Due to their low energy, secondary electrons are eventually released from the 
magnetic trap. Accelerated by the retarding potential, they hit the Si-PIN detector 
and thus produce background in the narrow energy interval of the signal /3-decay 
electrons [H]. 

Background bunches due to stored electrons from tritium beta decays were already 
observed in the Troitsk neutrino mass experiment [15j . Furthermore, measurements at 
the pre-spectrometer [16] and corresponding Monte Carlo (MC) simulations [17] revealed 
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that at a pressure level of p = 10~ n mbar a single nuclear decay can lead to an 
enhanced level of background of more than 3 ■ 10 -2 cps for several hours, which is well 
above the design limit of 10~ 2 cps. Of even greater importance is the fact that this 
background shows large non-Poissonian fluctuations, which can result in a significantly 
reduced neutrino mass sensitivity of KATRIN [H] . 

1.3. The ECR method 

To mitigate the background problems arising from stored electrons in the KATRIN 
spectrometers, a novel method to remove stored electrons must be implemented. This 
method should be capable of removing electrons with energies up to several hundred 
keV, it should work with a very small duty cycle, and it should not enhance other 
eventual sources of background. 

In the following we present a novel experimental method which has the promise 
to fulfill these requirements. It is based on stochastic heating of electrons by electron 
cyclotron resonance (ECR) [T8|[T9j. The ECR technique is based on the principle that 
the frequency /rf of an external high frequency (RF) field is adjusted to the cyclotron 
frequency f c of stored electrons. In our specific case the resonance region is set to the 
center part of the spectrometer. 

In case the resonance condition /rf = fc is met, the electron on average will gain 
a small amount of energy when passing the analyzing plane. Since a stored electron 
will pass the resonance zone 10 4 times within a 10 ms period of stochastic heating, its 
energy will be significantly increased, growing by more than one order of magnitude. 
Consequently, the ECR technique heats up electrons so that their cyclotron radius r c 
eventually becomes larger than the inner radius of the spectrometer walls (r s = 5 m), 
at which point they are absorbed. 

In the Mainz neutrino mass experiment it was shown that background from stored 
electrons could be reduced by applying high frequency pulses to an electrode of the 
spectrometer [20j[2T]. In this case, however, the electric field was mainly longitudinal to 
the magnetic field, and the frequency (1 — 2 MHz) was much smaller than the electron 
cyclotron frequency corresponding to the magnetic field at the electrode position, so 
that this background reduction method is significantly different from the ECR method 
presented in our paper. 

The paper is organized as following: first, we outline the basic principle of the ECR 
technique at a KATRIN spectrometer (section [2]). Then we discuss the experimental 
results obtained from test measurements of the ECR technique at the pre-spectrometer 
(section |3]), demonstrating its highly effective background reduction capabilities. Finally, 
we use extensive MC simulations to derive optimum operating parameters for the ECR 
method at the main spectrometer. 
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Figure 2. The top picture shows a 5 keV electron under stable storage conditions in 
the KATRIN main spectrometer over a time period of 2-10 -5 s. The electron trajectory 
is a superposition of the fast cyclotron motion, an axial oscillation, and a much slower 
magnetron drift. In the bottom picture an electron trajectory in the presence of an 
RF-field is shown. After about 1 ms, the electron's cyclotron radius becomes larger 
than the radius of the spectrometer and the electron is absorbed. The color coding 
denotes the kinetic energy (note that the energy scale in each picture is different). 
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2. The working principle of ECR at a KATRIN spectrometer 

The cyclotron frequency f c (B, 7) of a relativistic electron with Lorentz factor 7 is given 
by 

eB 



27TH17 

where B denotes the external magnetic guiding field, and e and m represent the absolute 
value of the electron charge and its mass. Sources of primary stored electrons include 
/3-decays of tritium molecules (HT, T2) and processes occurring during a-decays of 
the radon isotopes 219,220 Rn. The energy distribution of these primary electrons varies 
by more than five orders of magnitude from the eV level up to about one hundred 
keV [22TI28] . However, in terms of cyclotron frequency f c (B,'-f) this implies a variation 
of about 20% (1 < 7 < 1.2) only. The narrow frequency interval thus greatly facilitates 
the implementation of the RF-technique for removing electrons of very different energies. 

In the KATRIN spectrometers the magnetic field changes from its minimum value 
B c at the center to the ends of the vessel B e by four orders of magnitude (see figured]). 
Consequently, the f c (B) dependence easily compensates the small f c (l) variations, i.e. 
electrons of different energies will pass the resonance regions at different positions of the 
spectrometer. On the other hand, a magnetic field with steep gradients reduces the areas 
and times of resonance. Accordingly, we will show that by sweeping the RF-frequency 
the efficiency for stochastic heating of stored electrons can be improved significantly (see 
section |4~TT) . 

The net energy gain of an electron depends on the relative phase between the 
electron's cyclotron motion and the external RF-field. In the case of an RF-field with 
constant frequency applied in a constant magnetic field, no net energy gain can be 
achieved, due to a periodically shifting phase relation. At KATRIN, however, the 
phase between the electron cyclotron motion and the RF-field changes randomly: at 
each transition through the resonance-region the energy of the electron changes. As a 
consequence, the trapped electron will penetrate more or less deep into the high field 
regions, which in turn leads to a random phase change between the electron and the 
RF-field. 

In the KATRIN main spectrometer field layout, electrons with energies above about 
E max 100 keV cannot be trapped, as either their cyclotron radius at the analyzing 
plane becomes too large, or their storage is no longer stable due to non-adiabatic 
effects [13] . The stochastic heating of electrons thus has only to provide the energy 
difference between the starting energy and E max . 

Since the phase space for gaining energy is always larger than the constrained 
phase space for losing energy (the electron cannot have less than zero kinetic energy), 
a net energy gain is achieved after many resonance-passes. The basic principle of the 
ECR technique is illustrated in figure [2l It shows the trajectory of a 5 keV trapped 
electron for stable storage conditions and in the presence of an RF-field. The trajectories 
were simulated with the KATRIN simulation software for field calculation and particle 
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tracking Kassiopeia [HESIEQ]. 

3. Proof of principle at the KATRIN pre-spectrometer 

Based on the general considerations presented above, an experimental program was 
initiated to verify the background reduction capabilities of the ECR method at the 
KATRIN pre-spectrometer. In addition to a proof-of-principle demonstration, the 
measurements should establish that the background reduction can be achieved by rather 
short RF-pulses and that no other background sources are enhanced due to the RF-feed- 
in. On a more technical level, in view of the rather complex inner electrode systems of 
the spectrometers, the measurements should also demonstrate that the coupling of an 
external RF-field does not affect the integrity of these crucial system components. 

3.1. Experimental setup 

In the final KATRIN setup the pre-spectrometer affords the option to work as a pre- 
filter for /3-electrons, thereby reducing their flux by six orders of magnitude. In this case 
the pre-spectrometer acts as an electrostatic filter with a fixed retarding voltage set to a 
value of 300 V below the endpoint. At this operating point, the energy resolution of the 
pre-spectrometer of AE = 70 eV is sufficient to leave the interesting high-energy part 
of the /3-spectrum unmodified. Before being integrated in the full KATRIN setup, the 
pre-spectrometer was operated as a stand-alone facility with an extensive suite of test 
measurements [T^l[T6|[T7ll3 11433] . The major tasks of this program were the development 
of advanced technologies and experimental methods that later can be applied to the 
much larger main spectrometer. 

The pre-spectrometer (see figure [3]) has a length of 3.4 m and a diameter of 1.7 m. 
At both ends a superconducting magnet is installed providing a magnetic field of 4.5 T 
in the center of the solenoid and 15.6 mT in the center of the spectrometer. At one 
end an 8 x 8 Si-PIN diode array is mounted. As a novel design feature, when compared 
to the Mainz and Troitsk set-ups [20 |[34T436] . the vessel itself is set on (negative) high 
voltage. An inner electrode system consisting of two full electrodes and a wire electrode 
can be set to a different potential relative to the tank. 

The basic principle of the experimental setup for RF-injection at the pre-spectro- 
meter is shown in figure HI An RF-coupler box (see figure [5]) superimposes the RF- 
signal from an RF-power amplifier (model: TP 30/10, manufacturer: Schlumberger) to 
the DC voltages required to shape the electrostatic potential. In order to obtain high 
RF-voltages on the inner electrodes, it was examined whether resonant amplification 
could be obtained by exploiting the natural electric resonance of the inner electrode 
system. Although this electrode system had not been designed ab initio with respect to 
RF-properties, the presence of symmetric dipole-electrodes turns out to be favorable. A 
network analyzer was used to examine the impedance behavior of the electrode system. 
Low-impedance resonant frequencies were found at 32 MHz, 62.5 MHz and 89 MHz 
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Figure 3. Schematic view of the pre-spectromctcr. 1: superconducting solenoids. On 
the left end of the spectrometer an electron gun is installed to test the transmission 
properties [3T] (here the Pt30-2 krypton source is mounted); on the other end a 64- 
pixel Si-PIN diode detector is mounted, 2: pre-spectrometer vessel, 3: inner electrode 
system, 4: 90° pump port (housing the Au30-1 krypton source), 5: 45° pump port, 6: 
insulator 

(see figure E]). For the measurements reported on below, it was decided to use a fixed 
excitation frequency of /rf = 62.5 MHz. 

At this frequency, the dipole-electrodes act like a symmetric A/2 transmission line 
(see figure HJ). The pair of electrode connection rods form, by coincidence, a A/4 
transmission line feeding the A/2 dipole from one end, transforming the high impedance 
at this point into a low impedance at the feed point within the RF-coupler box. The 
narrow resonance peak of ~ 2 MHz width points to a high Q-factor, from which in turn 
a significant resonant amplification can be expected. The RF-power that was fed into 
the RF-coupler box was of the order of6-8W(17-20V RMS into 50 ft) only. 

With the RF-frequency being fixed at /rf = 62.5 MHz, the magnetic field at the 
center of the pre-spectrometer B^ s has to be adjusted to bring the electrons in resonance 
condition f c = /rf with the external RF-field. Figure [7] shows the resonance conditions 
for different magnetic field configurations (B^ s = 1,2,3 mT). The calculation shows 
that the largest resonance area is expected for B^ s = 2.15 mT. 

3.2. Simulation software 

Accompanying the measurements, MC simulations were performed using the simulation 
software Kassiopeia [291 [30], which has been developed over the past years by the 
KATRIN collaboration. Kassiopeia is a software package that comprises precise 
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Figure 4. Sketch of experimental setup: An RF-signal produced by an RF-gencrator 
is amplified and coupled into the inner electrode system of the pre-spectrometer, which 
can be operated in a dipolc mode, allowing the application of a transversal electric field. 




Figure 5. Electric circuit diagram of the RF-coupler box, which allows for a 
superposition of HV and RF-components to supply the inner electrodes. Left: 
Schematic view. Right: Photograph of the coupler box. 1) The "Balun" unit changes 
the single-ended signal to a differential signal (two opposite phase signals), it is 
constructed as a resonant circuit as shown in the schematics. 2) Capacitors to separate 
the HV DC component from the RF-input. 3) RF-choke and 4) capacitors working as 
LC low pass filter to prevent the RF-component from entering the HV supplies. 
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Figure 6. The measured signal reflection coefficient as a function of frequency. A 
dip in the reflection coefficient marks an area of resonance. A sharp resonance is an 
evidence of a high magnification of the electric field inside the spectrometer. At a 
frequency of 62.5 MHz (corresponding to a wavelength of 4.8 m) 70% of the applied 
energy was absorbed by the system. 



and fast calculation of electromagnetic fields and particle trajectories. The trajectory 
calculations of Kassiopeia are based on explicit Runge-Kutta methods described 
in [3TH39] . Electric and magnetic field calculations are performed via the zonal harmonic 
expansion [iQllJT] • In the case of electric fields, computations are carried out using the 
boundary element method |42j . 

Elastic, electronic excitation and ionization collisions of electrons with molecular 
hydrogen are included in the Kassiopeia. Total and differential cross section for elastic 
scattering are taken from [4~3] - l4"5] , data for inelastic scattering (electronic excitation and 
ionization) is based on [9HTT 4H&M^8] . Furthermore, arbitrary residual gas compositions 
of hydrogen, water, nitrogen or argon jl9l|50] can be chosen. The field, tracking and 
scattering simulations originate from FORTRAN and C codes developed between 2000 
and 2008 by one of us (F. G.). 

In the framework of the investigations presented below, Kassiopeia was equipped 
with a 83m Kr event generator taking into account the emission of conversion and 
Auger electrons, based on data in [5T],[52j. As a first approximation, the RF-field was 
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Figure 7. Cyclotron frequency of a 1 keV electron in the pre-spectrometer for three 
different values of B^ s . In the green colored regions the cyclotron frequency f c coincides 
with the frequency of the external RF-hcld /rf = 62.5 MHz. At a low magnetic field 
strength of B^ s = 1 mT this only occurs in narrow regions close to the entrance and 
exit of the spectrometer. When increasing B^ s to 2 mT, the area of resonance axially 
moves inward to the center and gets larger. Since electrons spend a large fraction 
of their axial oscillation period in the center and thus in resonance, a large effect is 
expected there. For magnetic field strengths of = 3 mT and higher, the area of 
resonance gets smaller and moves radially outwards before finally disappearing. 
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implemented in Kassiopeia as a sinusoidal function of the form 
E(x,t)=E (x) sin f fRF ^ " * \ , 

with Eq{x) and /rf denoting amplitude and frequency of the RF- field, respectively. The 
amplitude E (x) is chosen to be oriented perpendicular to the beam axis (E = \E \), 
since only this component of the field is able to increase the transversal energy E± of 
the electrons. The frequency can either be stationary /rf(£) = /rf, or swept through 
in steps /rf (t) = f^ cp (t)- 

The functionality and validity of the implemented RF-field routines in Kassiopeia 
were tested by comparing the numerical results to analytical calculations based on [19] 
and [18]. In this context, we consider parameters like the energy gain, position and 
momentum of low-energy electrons in a constant magnetic field under the influence of 
an RF-electric field. The agreement between all parameters produced by Kassiopeia 
and the analytic solution is of the order of 10~ 10 , which is of sufficient precision to 
validate the results presented in the following. 

3. 3. Source of stored electrons 

To investigate the efficiency of the ECR technique in background reduction, an intense 
source of high-energy stored electrons was required. In this work a gold-implanted 
(Au30-1) and a platinum-implanted (Pt30-2) rubidium-krypton ( 83 Rb/ 83m Kr) source 
were used (see figure [8]) [53]. The sources were installed at the horizontal pump port 
and at the electron gun position (see figure [3]). The implanted mother isotope rubidium 
decays into an excited state of krypton ( 83m Kr), which is metastable and has a lifetime 
of 1.83 h (see figure [8]). The emanation rate of 83m Kr of the Au30-1 source is c^ u = 11% 
and t£* = 6% for the Pt30-2 source [54] ■ 

The decay of 83m Kr produces conversion electrons in the multi-keV range (see decay 
scheme in figure [8j), leaving a vacancy in the electron shell. Subsequent relaxation 
processes lead to the emission of Auger electrons in an energy range of 10 2 — 10 4 eV. 
As mentioned above, the creation of conversion electrons as well as full subsequent 
relaxation processes involving the emission of several Auger electrons in one decay are 
included in Kassiopeia. The average multiplicity per decay is 11.6 ± 2.5 electrons in 
our model. 

The number of 83m Kr decays per second A-^r of both sources in the sensitive 
spectrometer volume was determined by a comparison of the measurement of the 
background rate at different central magnetic fields B^ s with a corresponding simulation 
(see figure [9]). In the latter, the number of electrons reaching the detector per 83m Kr 
decay was simulated for different values of B^ s . Since Akt is independent of B^ , the 
simulated relative rate was fitted to the measured rate with A^ T as a free parameter. 
With a x 2 minimization the best fit value was found at A^ r = 0.13 ± 0.002 s^ 1 . 

The resulting background level of &Kr = 6 cps at B^ s = 2.15 mT is significantly 
larger than the normal background rate of 6 intr ~ 10~ 2 cps without the krypton source. 
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Figure 8. Solid 83 Rb/ 83m Kr source. Left: Photograph of the platinum-implanted 
83 Rb/ 83m Kr source. Right: Decay scheme of 83 Rb. 



The ECR technique is thus targeted at removing these stored electrons by stochastically 
heating them up. 

3-4- Measurement results 

In the following we describe the experimental results obtained at the pre-spectrometer 
operated in a mode where an external RF-field is applied to the inner electrode system. 
The objective of the first measurement was to verify a strong reduction of the background 
rate when adjusting the magnetic field B^ to the resonance condition f c = /pp. The 
second measurement was targeted at showing that the ECR technique is indeed removing 
the primary stored electrons. The final goal was to prove that the ECR technique does 
not amplify other potential sources of background. 

3.4-1 ■ Resonance effect The objective of this measurement was to demonstrate a 
resonance behavior of the background reduction R centered at f c = / RF , with R defined 

as 



where Rrf, off denotes the background rate when no RF-field is applied and Rrf, on 
represents the rate with the RF-field applied. To find the maximal R, the magnetic field 
B^ s and thereby the cyclotron frequency f c were increased from a small value up to the 
expected resonance at B res = 2.15 mT and above. 

Figure [10] shows the experimental results of the background reduction factors R. 
The data feature a clear resonance centered at B TCS = 2.15 mT, as expected for ECR. 
The asymmetric shape of the resonance curve can be explained by the inhomogeneity 
of the magnetic field (see figure [7]). Its minimum value B^ s is reached in the center of 
the spectrometer, while increasing towards both ends. Correspondingly, if the central 
magnetic field is too low for resonance {B^ s < B TeB ), smaller resonant regions will exist 
only further away from the center. If the magnetic field is increased, these regions will 
move inwards and get larger, becoming maximal at the resonance field. However, if the 



R 



RrF, off — RrF, 
RrF, on 



on 



(2) 



3 PROOF OF PRINCIPLE AT THE KATRIN PRE-SPECTROMETER 



14 



77 60 

Q. 

U 

B 

£50 — 

"O 

rs 
o 



o 

CO 

CD 



30 — 



20 — 



10 — 



| — Measurement 

Kassiopeia simulation 



Ql I ± I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

2 4 6 8 10 12 14 16 

Magnetic field B PS [mT] 



CO 



10 §- 



10* 



10 H 



stored electrons at 2 mT 
stored electrons at 1 5 mT 




5000 10000 15000 20000 25000 30000 35000 

Energy[eV] 



Figure 9. Top: Measurement (full dots) and simulation (open dots) of the background 
rate for different central magnetic field strengths B^ s . The drastic drop of the rate with 
decreasing magnetic field is due to non-adiabatic effects. They occur when the magnetic 
field changes significantly within one gyration. In this case, the transformation of 
E± into E\\ and vice versa is no longer proportional to the change of the magnetic 
field. Consequently, the electron's polar angle changes randomly and eventually hits a 
value below the trapping threshold. Bottom: Energy spectrum of stored electrons at 
j3,? s = 15 mT compared to B^ s = 2 mT. From this plot it is evident that high energy 
electrons are affected most strongly by non-adiabatic effects: At low magnetic field, 
the number of high energy stored electrons is reduced by a factor of 10, whereas all 
low energy electrons are still stored. 
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central magnetic field is too high above the resonance value in the center (B^ s 3> B Tes ), 
the resonance condition f c = f^p is not met in the entire spectrometer volume. This 
explains the asymmetric shape of the curve above resonance. The finite width of the 
cut-off above B res results from the different energies of stored particles. The scan of -B^ s 
through the resonance was performed with two different RF- amplitudes Eq. In both 
cases an identical resonance behavior was observed (see figure [TO]) . 

The reduction factors as well as the characteristic resonance pattern for low values 
of B^ s and the asymmetric shape are reproduced by corresponding MC simulations 
using Kassiopeia [29l|30]. These simulations were also used to deduce the amplitude 
Eq of the RF-field in the pre-spectrometer, which could not be determined directly. 
As visible in figure HOj an amplitude of Eq = 100 V/m in the simulation yields 
comparable experimental reduction factors. This agreement implies actual field values 
of E = 100 - 200 V/m at full amplitude. 

Both measurement and MC simulations reveal that the fixed frequency setting used 
at the pre-spectrometer only removes part of the stored electrons (factors of 3 — 5), at 
resonance condition. The reason that not all stored electrons are removed is that we have 
used the technique of resonant amplification, which yields an RF-field of fixed frequency. 
Secondly, the rather steep gradients of the magnetic field of the pre-spectrometer are 
not optimal for achieving large resonant areas for electrons over a broad energy range. 
At the main spectrometer the magnetic field is more homogeneous over larger areas, 
thus we expect the ECR technique to be much more efficient there. 

3.4-2. Effect of RF -pulsing To test whether the ECR technique is indeed removing the 
primary stored electrons, we investigated the behaviour of the background rate after 
RF-pulsing. If the primary electrons are removed, we expect that after an RF-pulse 
the background rate returns back to its original value following an inverse exponential 
behavior R g = 1 — e . The parameter t T denotes the characteristic rise time of the 
rate, which is governed by the average storage time t stor of the trapped electrons and 
thus depends on the actual UHV conditions. 

Figure [TT] shows the measured behavior of R g as a function of the time t g since the 
last RF-pulse. An inverse exponential fit to the data yields t T = (9. 14 ±1.45) s, which is 
in good agreement with the simulated average storage time of tf™ T = (8.25±0.23) s. This 
value was achieved by a full MC simulation with Kassiopeia of 400 83m Kr decays in the 
pre-spectrometer at a pressure of p = 6 ■ 10~ 10 mbar with a realistic gas composition of 
50% H 2 , 25% H 2 and N 2 each. The parameter tf^ T was deduced from an exponential 
fit to the storage time distribution of all stored electrons. The error only includes the 
statistical error from the fit, whereas systematic uncertainties such as uncertainties in 
the pressure or the energy distribution are not considered here. 

As compared to the situation at the main spectrometer, typical storage times 
at the pre-spectrometer measurements are much shorter. This is due to the higher 
pressure, and secondly because mostly low energy electrons are stored in the magnetic 
field configuration of the pre-spectrometer with a central field of B^ = 2.15 mT (see 
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Figure 10. Measurement of the relative reduction factor R of the background from 
stored electrons as a function of the magnetic field strength in the center of the pre- 
spectrometer B^ s . As expected, for the resonance condition / c = /rf the maximal 
reduction is observed at B^ s = 2.15 mT. The two RF-field amplitudes Eq were 
controlled by adjusting the power fed into the system. The actual field inside the 
spectrometer can be determined by comparison to the simulation. The data show the 
expected resonance pattern of the ECR technique. 



figure |9]). The cross section for ionization is maximal for 10 2 eV electrons and decreases 
for higher energies |4Tj. Therefore the cooling for low energy electrons is faster. 

3.4-3. Influence on other background sources Apart from stored electrons from nuclear 
decays studied above, an important background mechanism is related to cosmic 
muons [32l [55] that hit the spectrometer vessel and generate low-energy electrons. 
Starting from the inner surface, these electrons are electrically and magnetically shielded 
from the sensitive volume of KATRIN [56], since the inner electrode is on a more negative 
potential than the vessel hull and the magnetic field lines run parallel to the electrode 
surface. To test whether an RF-field influences this important shielding effect, ECR 
measurements without the krypton sources were performed at = 2.15 mT. 

Without the krypton sources we measure an intrinsic background rate of bi ntv = 
(11.5±0.9)-10 -3 cps in the detector energy region of interest (ROI: 15 keV - 21 keV). The 
background rate was reduced to &i ntr = (6.2 ±1.2) ■ 10" 3 cps when applying a permanent 
RF-field. Since the pre-spectrometer background rate largely originates from the a- 
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Figure 11. Measurement of the increase of rate after 3 s long RF-pulses at a pressure 
of p = (6 ± 2) ■ 10~ 10 mbar and a magnetic field of B^ s = 2.15 mT. Each data point 
represents a separate w 10 h measurement during which the RF-pulse was alternatingly 
switched on for i on = 3 s and switched off for t a g = 5—50 s, respectively. The plot shows 
the average rate during t a g. The rate R g{t s = s) of 1.37 ± 0.07 cps was measured 
with a continuous RF-ficld. The horizontal error bars indicate the inaccuracy of the 
timing, accordingly R g{t f[ = s) has no horizontal error bar. An inverse exponential 
fit yields t x = (9.14 ± 1.45) s. 



decays of 219 Rn, 220 Rn emanating from the NEG strips and structural materials [T61I32] . 
the reduction can likely be explained by the active removal of stored electrons following 
radon decays. 

From this result we can put a conservative upper limit on the cosmic muon induced 
background of 6 muon < (6.2 ± 1.2) • 10~ 3 cps at the pre-spectrometer, and we can exclude 
that this background is substantially increased by the application of the RF-field. This 
is an important result when discussing the use of the ECR method during the long-term 
scanning of the tritium /3-spectrum with the main spectrometer. 



4. Expected background reduction at the main spectrometer 

To validate and study in detail the concept of stochastic heating by the ECR method as 
a means of background reduction at the main spectrometer, extensive MC simulations 
were performed. These investigations are motivated by the expected background rates 
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from stored electrons obtained by extrapolating the experimental results observed at 
the pre-spectrometer [16] to the significantly larger main spectrometer (diameter = 
10 m, length = 24 m) [TJ]. In the following we present the results of MC simulations 
with Kassiopeia that probe the efficiency of the ECR technique in alleviating this 
background at the main spectrometer. 

Kassiopeia was used to implement the full electromagnetic design of the main 
spectrometer, consisting of a system of superconducting coils and an external air coil 
system with a diameter of d = 12 m, which together provide an adiabatic guidance 
system for signal /3-decay electrons through the spectrometers. The precision retarding 
potential is created in the following way: the spectrometer vessel itself is on negative 
high voltage (HV) and a wire-based inner electrode system is operated on a slightly 
more negative HV, allowing for a 10 -6 precision of the filter potential [56]. The RF-field 
was implemented as described in subsection I3.2L 

The aims of the simulations described below are to optimize the frequency settings, 
to test the efficiency of ECR with low field amplitudes, to investigate the effects of 
ECR on electrons of different energies, and finally to analyze the impact of ECR on the 
statistical uncertainties in the neutrino mass analysis of KATRIN. 

4-1- Optimizing the frequency settings 

To compare the two possible frequency settings and f^ cp ) as well as to optimize 
the sweeping strategy we first carried out simulations with a simplified configuration 
in a constant magnetic field -B C onst = 0.3 mT (matching the magnetic field in the 
center of the main spectrometer). The starting energy of the electrons was fixed to 
-Kprim = 5 keV (matching the typical energies of electrons from nuclear decays). The 
cyclotron frequency of these electrons is f c = 6 MHz, which corresponds to a cyclotron 
period of t c = 1.67 • 10 -7 s. 

As expected, the simulation demonstrates that no net energy gain can be achieved 
with a constant frequency RF-field in a constant magnetic field. Instead, the energy 
of the electron changes periodically with time (see figure [T2I) . In a constant magnetic 
field, a sweeping frequency mode is needed to heat up electrons. The frequency sweep 
is implemented as a stepwise increasing function, where the duration of a step t stC p and 
the number of steps n stcp within a sweep are free parameters. The simulations visualized 
in figure [12] show that the energy change of the electron by ECR is most efficient, if 

• the duration of a single step is of the order of a few times the cyclotron period, i.e. 

^step ^ Tl ' t c . 

• n step is large so that many different frequencies occur within one sweep. 

Following these initial results for a fixed energy and a constant magnetic field, we 
expand our investigations to the realistic case of the main spectrometer. The magnetic 
field there increases axially towards the field-generating solenoids. Consequently, the 
cyclotron frequency f c of trapped electrons changes along their trajectory. Therefore, 
an RF-field at a fixed frequency can lead to a net energy gain. By varying the frequency 
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Figure 12. Comparison of example frequency settings. The figure shows the kinetic 
energy of an electron with a starting energy E pi i m = 5 keV in a constant magnetic 
field of -Bconst = 0.3 mT in the presence of an RF-field with constant or sweeping 
frequency mode and an amplitude of E = 100 V/m as a function of time. The 
constant frequency was set to / CO nst = 8.28 MHz. This case (green) leads to a periodic 
energy change. In the sweeping mode the frequency was varied between / m i n = 7 MHz 
and / m ax = 8.5 MHz. The frequency is swept through in n stcp steps with an individual 
time step £ s t e p- If the time steps are too short i s t ep = 0.01 fis (violet, +20 keV) or 
the number of sweeping steps is too small n s t op = 15 (red, +40 keV), no significant 
energy increase can be achieved. The best result is obtained when the scanning of the 
frequency occurs in small steps with a duration of a few times the cyclotron period 
(blue, +60 keV). The energy curves for sweeping RF-cases have been separated by 
offsets of 20 keV, 40 keV and 60 keV, respectively, for better visibility. 

we expect to further improve the background reduction efficiency. This was investigated 
by considering the test case of electrons with energies up to 18.6 keV from tritium (3- 
decay in the main spectrometer volume, both with constant and sweeping frequency of 
the RF-field. 

The constant RF-frequency is chosen to correspond to the cyclotron frequency 
f c of a low-energy electron (1 eV) in the minimum magnetic field B^ s . In the sweeping 
frequency mode, the frequency f^ ep is changed between 8 and 10 MHz, corresponding 
to cyclotron frequencies of 1 eV and 100 keV electrons at B^f , respectively. The 
parameters of the RF-field settings are summarized in table [TJ 

The comparison of the two RF-modes shows that in case of a sweeping frequency, all 
electrons could be removed, whereas in the case of constant frequency a fraction of 6.5% 
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Table 1. RF- field settings for the simulation with Kassiopeia 
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Table 2. Percentage of electrons remaining trapped in the presence of an RF- field with 
different settings. The cases of no RF-fields and RF-fields with constant or sweeping 
frequency are compared. For these investigations an ensemble of 2 • 10 3 tritium /3-decay 
electrons was simulated for each RF-configuration. Electrons are considered to be in 
stable storage condition if their storage time exceeds 500 ms. Without the RF-field, 
a fraction of 51.5% of the /3-electrons are stored, the remainders immediately hit the 
wall due to their large cyclotron radius, or leave the spectrometer axially due to non- 
adiabatic effects [Tl]. The table demonstrates the high efficiency of the ECR technique, 
in particular in the sweeping mode with a complete reduction of stored electrons. 



RF-setting 


Remaining percentage of trapped electrons 


off 


51.5% 


constant frequency 


6.5% 


sweeping frequency 


0% 



of the electrons remains trapped within a fixed time interval of 500 ms (see table [2]). 
In case that not all electrons are removed within a single RF-period, the subsequent 
RF-pulses however will remove a significant fraction of these trapped electrons. The 
background level in between subsequent RF-pulses is thus unaffected by these electrons, 
which remain "parked" at high energies, as the energy- dependent cross section for 
ionizing collisions decreases in the keV-range |47j . 

4-2. Optimizing the amplitude of the RF-field 

In order to comply with the long-term integrity of the main spectrometer wire electrode 
system with its 24000 wires of a diameter of 300 /im (outer layer) and 200 fim (inner 
layer) [56], the lowest possible amplitude of the RF-field E is mandatory. For low E 
values, the electrons have to pass the analyzing plane more often to reach the extraction 
energy of E max « 100 keV and, consequently, longer removal times are expected for 
lower amplitudes. 

To test this hypothesis for trapped electrons from tritium /3-decay, the effect of a 
sweeping RF-field with 20 different amplitudes between Eq = 10 — 200 V/m was tested 
(the settings of the sweeping frequency are given in tabled]). As an important result 
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Figure 13. Removal times t lem of tritium /3-decay electrons in the presence of an 
RF- field with sweeping frequency as a function of RF-amplitude Eq. For this figure, 
10 tritium /3-dccay electrons were simulated for each Eq. The parameter t rem is given 
by the lifetime r, which was determined by fitting an exponential curve to the storage 
time distribution of the electrons for each Eq . The figure shows that i rem decreases for 
larger values of Eq. However, even for very low values of Eq = 10 V/m the removal 
time i rem remains smaller than 10 ms. 



we find that for the values of E considered all stored electrons could be removed. The 
average removal time increases by about an order of magnitude for smaller values of Eq, 
as expected (see figure [T3|) . However, even for very small amplitudes of Eq = 10 V/m 
the removal time is still smaller than 10 ms. 



4-3. Effect of an RF-field on electrons of different energies 

During the long-term tritium scanning measurements the background of stored electrons 
will not only result from the above considered tritium /3-decays in the spectrometer 
volume, but will also comprise electrons from 219,220 Rn a-decays with energies of up 
to 450 keV [23l[2l]. Consequently, we investigate the effectiveness of an RF-field as a 
function of the initial electron energy E wmi . 

For this investigation, electrons of nine different starting energies in a range between 
10 eV — 100 keV were simulated in the presence of an RF-field with sweeping 



E 



prim 
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frequency (the parameters are given in tabled]). For each of these values, an ensemble 
of 2 • 10 4 electrons was started isotropically in the main spectrometer. Again, the ECR 
method works efficiently, as even for this extended energy range all stored electrons were 
removed within 1 ms. This demonstrates that the ECR method is a promising tool to 
remove both electrons from tritium /3-decays as well as from 219 ' 220 Rn a-decays. 

4-4- Impact of ECR on statistical uncertainty of KATRIN 

In [H] it was demonstrated that the background due to stored electrons following 
219 ' 220 Rn a-decays and tritium /3-decays results in a significant increase of the statistical 
uncertainty cx"^ 1 of the observable m| of KATRIN. In the following we investigate how 
efficiently the application of the ECR technique can reduce af^ and thereby improve 
the statistical neutrino mass sensitivity. 

As shown in the sections above, the ECR technique can be characterized by an 
almost negligible removal time t rem < 10 ms and a 100% efficiency of removing stored 
electrons of all energies. The only free parameter that influences the background level 
from nuclear decays is thus the time interval between two subsequent RF-pulses At pu i se . 
Since during the RF-pulsing no neutrino mass data can be taken, one aims to maximize 
Atp U i sc . Accordingly, we present the background level b nuc \ arising from nuclear decays 
and cr^ 1 as a function of At pu i se (see figure [T4|) . 

Due to the good UHV conditions at the main spectrometer (p = 10~ n mbar), 
the ionization times t- lon of keV electrons are rather long (e.g. for a 5 keV electron 
^ion ~ 3 min at p — 10~ n mbar). We thus expect that correspondingly large values of 
Atp U i se are sufficient to reduce substantially the background rate. 

For the statistical analysis we use the detailed background model described in [T3] . 
which is based on MC simulations of 219 ' 220 Rn a-decays and tritium /3-decays in the 
main spectrometer and experimental data in [16]. We expect a background rate from 
these decays of b QUC \ = 3 ■ 10~ 2 cps [H] . 

The statistical uncertainty cr^ L is determined by a \ 2 fit of the theoretical integral 
/3-spectrum for 5 • 10 3 simulated KATRIN measurements, considering the nuclear decay 
model as the only background source. For each simulation m Ve = eV is assumed. The 
width of the distribution of the fitted neutrino mass squared nip determines crf^'- The 
statistical analysis, visualized in figure [HI reveals following facts: 

• With a value of At pu i sc = 2 min the background from stored electrons can be 
eliminated almost completely to a level of 6 nuc i = 6 • 10 -4 cps. The statistical 
uncertainty of cr^"f' = 0.011 associated with this background is comparable to 
the situation without background where the statistical error is dominated by the 
variation of the signal count rate. With a pulse length of t p = 10 ms and a repetition 
time of Atp U i se = 2 min a duty cycle of t p j (At pu \ sc + t p ) ~ 8 • 10~ 5 would be reached, 
which implies a negligible loss of neutrino mass measurement time. Moreover, this 
RF-pulsing rate ideally fits to the measurement schedule of tritium scanning [1] 
where the retarding potential is varied on the order of minutes. Accordingly, the 
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Figure 14. The figure shows the background arising from nuclear decays 6 nuc i 
(left axis, circles) and the associated statistical uncertainty er" t " c t h of the observable 
m§ (right axis, triangles) as a function of the time between RF-pulses Ai pu i sc . A 
background model of 3 • 10~ 2 cps without ECR is assumed. A value of At pu i so = 2 min 
is sufficient to reduce the background to less than 10 -3 cps. 



RF-pulses could be applied while setting the retarding potential. 

• Repetition times of At pu i sc = 10 min reduce the original background level of 
&nuci = 3- 10~ 2 cps to 6 nuc i = 4.3- 10~ 3 cps, corresponding to a statistical uncertainty 
of cr^t 1 = 0.016. As a comparison, this uncertainty is of the same size as the 
statistical uncertainty associated with a constant (Poisson distributed) background 
level of 10 -2 cps as anticipated in [I]. 

• For long repetition times At pu i se > 6 h the ECR technique is no longer effective. 

These considerations are also valid for all other active background reduction mechanisms 
that allow for complete removal of all trapped particles in time periods of the order of 
ms. 

5. Conclusion 

Previous measurements [16] and simulations [HJ[T7] have revealed that stored electrons 
arising from nuclear decays in the volume of the KATRIN spectrometers can cause 
large background rates that exceed the design background limit of 10 -2 cps. It was 
shown that the non-Poissonian nature of this background can significantly influence the 
neutrino mass sensitivity of KATRIN. 

In this paper, we have demonstrated the feasibility of a novel background reduction 
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technique based on stochastic heating of stored electrons via the process of electron 
cyclotron resonance (ECR). 

Measurements at the KATRIN pre-spectrometer have proven the basic functionality 
of the ECR method. In particular, we have demonstrated that rather short RF-pulses 
of At < 2 s are sufficient to achieve a large background reduction factor of about five. 
Moreover, we could show that the RF-field does not increase backgrounds from other 
sources, in particular muon-induced background. 

In addition to these measurements, we have carried out extensive simulations of 
the ECR technique with the Kassiopeia code to study the efficiency of background 
reduction at the main spectrometer. Two RF-operation modes were considered: a mode 
with constant frequency, as well as a mode with sweeping frequency. The simulations 
revealed the sweeping frequency mode yields a complete background reduction, while 
the technologically simpler constant frequency mode still removes a large part (93.5%) of 
stored electrons. Moreover, the simulations with Kassiopeia proved that the stochastic 
heating is efficient for electrons with energies up to 100 keV, thus covering the entire 
energy range of stored electrons from tritium /3-decays and 219,220 Rn a-decays. 

Most importantly, the simulations showed that even in the case of a very low RF- 
amplitude of only E = 10 V/m, which is preferable with respect to the integrity 
of the wire electrode, all stored electrons can be removed within less than 10 ms. 
Furthermore, statistical simulations revealed that the time between two subsequent RF- 
pulses of Atp U i se « 2 min suffices to almost completely eliminate the contribution from 
the background due to stored electrons to the statistical error budget of KATRIN. The 
low RF- amplitudes and the small duty cycle of ~ 8 • 10~ 5 are important prerequisites 
to assure the integrity of the wire electrode system while at the same time not affecting 
the measuring time of KATRIN. 

Based on the experimental work at the pre-spectrometer presented here, as well as 
on the corresponding MC simulations described above, we intend to test and optimize the 
ECR method, as well as other active reduction methods at the upcoming commissioning 
measurements at the KATRIN main spectrometer, with the aim of operating the 
spectrometer nearly free of background. 
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